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ABSTRACT 

In this paper, we report on our analysis using Hubble Space Telescope astrometry and Kcck-I 
HIRES spectroscopy of the central six stars of Tycho's supernova remnant (SN 1572). With 
these data, we measured the proper motions, radial velocities, rotational velocities, and chemical 
abundances of these objects. Regarding the chemical abundances, wc do not confirm the unusu- 
ally high [Ni/Fe] ratio previously reported for Tycho-G. Rather, we find that for all metrics in all 
stars, none exhibit the characteristics expected from traditional SN la single-degenerate-scenario 
calculations. The only possible exception is Tycho-B, a rare, metal-poor A-type star; however, we 
are unable to find a suitable scenario for it. Thus, we suggest that SN 1572 cannot be explained 
by the standard single-degenerate model. 

1. Introduction 
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Type la supcrnovae (SNe la) are of great in- 
terest. They represent some of the most extreme 
physical situations in stellar astronomy, control 
the chemical evolution of galaxies and the Universe 
at intermediate to late times by producing large 
amounts of iron-group elements, and are uniquely 
powerful cosmic distance probes. But despite their 
wide-ranging significance, fundamental uncertain- 
ties remain around the progenitors of these cata- 
clysmic events. 

There is general consensus that SNe la are 
caused by the deflagration/detonation of a car- 
bon/oxygen white dwarf (CO- WD) which is ac- 
creting material from a binary companion. Sce- 
narios exist where the explosion can be initiated 
from a detonation on the surface of the star (Livne 
& Arnett 1995; Fink et al. 2010), through runaway 
carbon burning in the white dwarf's interior, or 
through a cataclysmic merger of objects. 

Observationally, two main models for this ac- 
cretion process can be identified. The single- 
degenerate scenario (SD-scenario) sees the accre- 
tion process occurring through Roche-lobe over- 
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flow (RLOF) of a close nondcgcncrate compan- 
ion (also known as a donor star). This compan- 
ion, which has undergone common-envelope evolu- 
tion with the white dwarf, can be a helium, main- 
sequence, subgiant, or red giant star. In all cases 
the donor star should survive the explosion (ex- 
cept for possibly in the case of the helium-star 
donor; R. Pakmor 2012, private communication) 
and remain visible post-explosion. 

The second scenario is the dynamical merger of 
two white dwarfs, known as the double-degenerate 
scenario (DD-scenario). In this scenario, the co- 
evolution of two stars eventually leads to a close bi- 
nary of two white dwarfs, which are able, through 
the emission of gravitational radiation, to merge 
over a wide range of times after the initial forma- 
tion of the system. In most cases this would leave 
no remaining star (e.g., Pakmor et al. 2010). 

Both scenarios have support in observations 
and theory. The detection of circumstellar ma- 
terial around certain SNe la (Patat et al. 2007; Si- 
mon et al. 2009; Sternberg et al. 2011; Foley et al. 
2012) provides evidence for the SD-Scenario. On 
the other hand, the lack of substantial hydrogen in 
the majority of other SNe la (Leonard 2007) poses 
a challenge to the SD-Scenario. 

Kasen (2010) suggests that the interaction of 
the shock wave with the nondegenerate compan- 
ion should result in a light excess at early times of 
an SN la light curve, which depends on the viewing 
angle and the companion radius. Such an excess 
has not yet been observed (Hayden et al. 2010; 
Tucker 2011; Bianco et al. 2011; Ganeshalingam 
et al. 2011), which is at odds with red giant com- 
panions forming the majority of SNe la. Justham 
(2011), Di Stefano et al. (2011), and Hachisu 
et al. (2012a,b), however, have suggested a sce- 
nario where the white dwarf is spinning and thus 
can accrete above the Chandrasekhar limit. The 
explosion would only occur once the white dwarf 
had spun down sufficiently, which would give the 
red giant a chance to evolve and would not require 
the detection of the early excess in the light curve 
in a red giant progenitor scenario. 

Population-synthesis calculations are challeng- 
ing, with various authors getting different results 
for the same inputs. However, there is a general 
trend from these calculations that neither single- 
degenerate nor double-degenerate stars can pro- 
vide enough systems to explain the observed SN la 



rate (Han 2008; Ruiter et al. 2009; Mennekens 
et al. 2010; Yu & Jeffery 2010). Several authors 
suggest that the population might comprise both 
single-degenerate and double-degenerate systems. 

The physics of white-dwarf mergers is challeng- 
ing to simulate numerically, but in the simplest 
calculations, these mergers lead to the formation 
of a neutron star via electron capture, rather than 
to a thermonuclear explosion (Saio & Nomoto 
1985). Recently, Pakmor et al. (2010) have shown 
that for certain parameters (white-dwarf binaries 
with a mass ratio very close to 1) the merger may 
explain subluminous supernovae (e.g., SN 1991bg; 
see Filippenko 1997, for a review), although Dan 
et al. (2011) note that the initial conditions of the 
system may change these conclusions. 

SN 2011fe was detected only ^ 11 hr after the 
explosion, and (with a distance of 6.4 Mpc) is one 
of the closest SNe la ever found (Nugent et al. 
2011). Nugent et al. (2011) and Brown et al. 
(2011) have not found any early-time light-curve 
excess predicted by Kasen (2010) and thus rule out 
a red giant donor. Radio and X-ray observations 
by Horesh et al. (2011) show no strong signs of 
pre-explosion outflows, which again contradicts a 
red giant scenario for SN 201 Ife. Additional radio 
measurements by Chomiuk et al. (2012) suggest 
a low density around SN 201 Ife, which is at odds 
with many conventional single-degenerate scenar- 
ios. Li et al. (2011) have searched pre-explosion 
archival images and can also rule out luminous red 
giants and almost all helium stars as donors. In 
addition, Bloom et al. (2012) have used images 
believed to have been taken 4hr post-explosion 
and suggest the companion radius to be less than 
0.1 R0. Most of these results are consistent with 
a main-sequence companion or white-dwarf com- 
panion. 

Because it is very difficult to obtain robust con- 
straints on the progenitor system in the immedi- 
ate aftermath of a 10''''^ erg explosion, an alterna- 
tive is to study somewhat older and more nearby 
SNe that can be observed in great detail. Ruiz- 
Lapuente et al. (2004, henceforth RP04) have tried 
to directly detect donor stars in old and nearby 
SN la remnants within the Milky Way. They have 
identified two historical Galactic SNe well suited 
to this task — SN 1006 and SN 1572 (Tycho's 
SN). Both remnants are young (1000 and 440 yr 
old, respectively), closeby (2.2 ± 0.08 kpc, Winkler 
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et al. 2003; 2.8 ± 0.8 kpc, Ruiz-Lapucntc 2004, re- 
spectively), almost certainly SNe la from their ob- 
servational signatures (Badenes et al. 2006; Ruiz- 
Lapuente 2004; Krausc ct al. 2008; Rest et al. 
2008), and not overwhelmed by Galactic extinc- 
tion. In this paper, we will focus on SN 1572. 

RP04 investigated most bright stars in the cen- 
tral region of SN 1572 and found a star with an 
unusual spatial motion (Tycho-G by their nomen- 
clature); they suggested this as a possible donor 
star for SN 1572. While the star has an imusual 
spatial motion compared to other stars in the field, 
its current location and proper motion place it a 
significant distance from the center of the super- 
nova remnant (SNR) — a feature difficult to ex- 
plain in connecting Tycho-G to SNR 1572. 

One consequence of RLOF is a rotational ve- 
locity induced on the donor star by tidal lock- 
ing in the system. This results in an unusually 
large rotational velocity, related to the orbital ve- 
locity of the binary system, and it can be used 
to single out possible donor stars from nearby un- 
related stars. Kerzendorf et al. (2009, henceforth 
WEK09) investigated the rotation of Tycho-G but 
found no excess rotational velocity compared to a 
normal star. A comparison of WEK09's measure- 
ments of Tycho-G, including a revised radial veloc- 
ity frad, with Galactic kinematic models showed 
that it is statistically consistent with an unrelated 
thick/thin-disk star. However, WEK09 were able 
to provide an a priori unlikely donor-star scenario, 
where the star was able to lose its rotational sig- 
nature. 

Gonzalez Hernandez et al. (2009, henceforth 
GH09) analyzed a spectrum of Tycho-G observed 
with the High Resolution Echelle Spectrometer 
(HIRES; Vogt et al. 1994) instrument on the Keck- 
I 10 m telescope, finding a Wrad value consistent 
with WEK09's revised v^ad- They also measured 
the stellar parameters and metallicity of Tycho-G, 
concluding that it has an unusually high nickel 
abundance, which they claim can be attributed to 
the accretion of ejecta material on the donor star 
during the explosion. 

In this paper we analyze HIRES spectra of 
the six bright stars near the center of SNR 1572. 
These spectra were taken as part of the same pro- 
gram that obtained the data used by GH09, and 
we independently reanalyze the Tycho-G spec- 
trum in our program. We describe the observa- 



tional data and our data-reduction procedures in 
§2. Section 3 is divided into six subsections de- 
tailing the measurements of proper motion, ra- 
dial velocity, rotation, stellar parameters, and 
abundances, and we provide a detailed compar- 
ison between our and GH09's measurements of 
Tycho-G. In §4 we analyze the measurements of 
each star to investigate its potential association 
with SNR 1572, and we present our conclusions in 
§5. 

2. Observations and Data Reduction 

We obtained spectra with the HIRES spectro- 
graph on the Keck-I telescope on Mauna Kea. The 
observations were made on 2006 September 10 
and 2006 October 11 UT. Slits B5 and CI (width 
0.86"; B5 length 3.5", CI length 7.0") were used, 
resulting in wavelength coverage of 3930-5330 A, 
5380-6920 A, and 6980-8560 A with R = A/AA w 
50,000, providing us with the necessary spectral 
resolution and wavelength coverage to determine 
stellar parameters. 

The spectra were reduced using the makee 
package. All spectra were corrected to heliocen- 
tric velocities, using the makee sky-line method. 
The spectra were not corrected for telluric absorp- 
tion lines, but only regions known to be free from 
telluric contamination were used in the analysis to 
derive the stellar parameters. The final exposure 
times of the combined spectra for each candidate 
and the signal-to-noise ratio (S/N) at 5800-5900 A 
are shown in Table 1. Finally, we normalized the 
spectrum using the IRAF^ task continuum. We 
note that Tycho-C and Tycho-D were observed on 
the same slit (CI); they are separated by 2.1", 
and the seeing was ^ 0.8", with Tycho-C being 
roughly five times brighter than Tycho-D. 

In addition, we obtained low-resolution spec- 
tra {R ss 1200) of Tycho-B with the dual- arm 
Low- Resolution Imaging Spectrometer (LRIS; Oke 
et al. 1995) mounted on the Keck-I telescope. The 
data were taken on 2010 November 7 UT, using 
only the blue arm with the 600/4000 grism and 
the l"-wide slit. This resulted in a wavelength 



IRAF; the Image Reduction and Analysis Facility is dis- 
tributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Re- 
search in Astronomy (AURA) under cooperative agreement 
with the National Science Foundation (NSF). 
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Table 1 
Observations 



Tycho 


a(J2000) 


J(J2000) 


Date 


Slit 


texp V- S/Nb 


(Name) 


(hh:min:ss.ss) 


(dd:mm:ss.ss) 


(dd/mm/yy) 




(s) (mag) 



A 


00:25:19.73 


+64:08:19.60 


10/09/06 


B5 


900 


13,29 


~ 48 


B 


00:25:19.95 


+64:08:17.11 


10/09/06 


B5 


1200 


15.41 


~45 


C 


00:25:20.40 


+64:08:12.32 


11/10/06 


CI 


10,800 


19.06'= 


- 8 


D 


00:25:20.60 


+64:08:10.82 


11/10/06 


CI 


10,800 


20.70 


~ 3 


E 


00:25:18.29 


+64:08:16.12 


11/10/06 


CI 


9000 


19.79 


~ 9 


G 


00:25:23.58 


+64:08:02.06 


10/09/06 & 11/10/06 


B5&C1 


24,000 


18.71 


~ 25 



'^Magnitudes from RP04. 

''The S /N value was obtained by measuring the root-mean square of the pixels (each resolution element 
is sampled by 2 pixels) in continuum regions near 5800-5900A. For the purposes of measuring the stellar 
parameters, the spectrum was convolved and the S/N increases by a factor of 2.24. 

'=RP04 notes that this is an unresolved pair, with a brighter bluer component {V = 19.38 mag) and a 
fainter redder component {V = 20.53 mag). 



coverage of 3200 5600 A. These data were taken 
to obtain a precise measurement of the surface 
gravity for Tycho-B using the size of the Balmer 
decrement (Bessell 2007). 

The spectrum of Tycho-B was reduced using 
standard techniques (e.g., Foley et al. 2003). Rou- 
tine CCD processing and spectrum extraction 
were completed with iraf, and the data were 
extracted with the optimal algorithm of Horne 
(1986). Wc; obtaincxl the; wavelength scale from 
low-order polynomial fits to the calibration-lamp 
spectra. Small wavelength shifts were then ap- 
plied to the data after measuring the offset by 
cross-correlating a template sky to the night-sky 
lines that were extracted with the star. Using 
our own IDL routines, we fit a spectrophotometric 
standard-star spectrum to the data in order to 
flux calibrate Tycho-B and remove telluric lines 
(Horne 1986; Matheson et al. 2000). 

3. Analysis 
3.1. Astrometry 

Proper motions can be used to identify poten- 
tial donor stars because donor stars freely travel 
with their orbital velocity after the SN explosion 



disrupts the system. RP04 suggested Tycho-G as 
a possible donor due to its unusually high val- 
ues for both the proper motion and radial ve- 
locity. For this work wc measured proper mo- 
tions for 201 stars within one arcminute of the 
remnant's center. We used archival Hubble Space 
Telescope (HST) images for three different epochs 
(iJST Programs GO-9729 and GO-10098; Novem- 
ber 2003, August 2004, May 2005), each con- 
sisting of three exposures (Is, 30s, and 1440s) 
with the F555W filter using the Advanced Cam- 
era for Surveys (ACS). The scale in each exposure 
is 50 mas pixel" ^. This dataset results in a maxi- 
mum baseline of 18 months. 

We used an image from the middle epoch (2004) 
to establish a reference frame and oriented the 
pixel coordinate system with the equatorial sys- 
tem. We then applied a distortion correction for 
the F555W filter (Anderson & King 2006) and cal- 
culated transformations between all other images 
and the reference image. Next we used these trans- 
formations to calculate the position of all stars in 
the reference coordinate system, with the overall 
uncertainty of each position estimated. Some faint 
stars were not detected in the shorter exposures 
and were thus excluded from proper-motion mea- 
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surements. In total, 114 stars were used in the 
astrometric analysis. 

For each star, we fit a linear regression for the 
stellar positions over time in the pixel coordinates 
(which were aligned with the equatorial system). 
The X and y data were treated as independent 
measurements, with separate regressions solved 
for each axis direction. Uncertainties were esti- 
mated using standard least-squares analysis and 
the individual uncertainty estimates of each ob- 
ject's positions. 

There are three J2000 measurements of the ge- 
ometric center of SN 1572 from different datasets. 
Reynoso et al. (1997) used Very Large Array 
(VLA) data to measure the center as a = 00'*25™14?95, 
5 = +64°08'05.7"; Hughes (2000) used ROSAT 
data to measure a = 00^25'"19?G, 5 = +64°08'10"; 
and Warren et al. (2005) used Chandra data to 
measure a = 00^25™19f40, 5 = 64°08'13.98". We 
note that the X-ray centers agree rather well with 
a difference of less than 5", but the radio center is 
located roughly 30" away from the X-ray centers. 
Thus, we believe the error in the geometric center 
to be rather large (of order 30"). 

Table 2 lists the proper motions and uncer- 
tainties of all stars mentioned by RP04 (19 stars) 
which were analyzed in this work, as well as the 
distance to the geometric X-ray center measured 
by Chandra. 

We note that Tycho-2 has a relatively high 
proper motion, but its position in the year 1572 
was 67.95" away from the remnant's center, and 
we thus exclude it as a viable candidate for the 
donor. 

In Figure 1, we compare the distribution of 
proper motions of all measured stars to our candi- 
dates. All of our candidates are reconcilable with 
a normal proper-motion distribution. 

3.2. Radial Velocity 

For the radial- velocity measurement we first ob- 
tained well-calibrated wavelength solutions for our 
spectra. MAKEE performs an initial calibration 
of the wavelength using arcs and then refines it 
by cross-correlating the night-sky lines for each 
observation and determining minor offsets. Both 
science objects and radial-velocity standards were 
reduced in the same fashion. 

Each order of each star spectrum was then 



cross-correlated using the iraf task FXCOR (Tonry 
& Davis 1979) with at least two other radial- 
velocity standards (HR6349, HR6970, and HR1283) 
which had been observed on the same night. We 
measure the radial velocity of the standards and, 
comparing to the canonical values (Stefanik et al. 
1999), wc obtain a systematic error of ^ 1 km s~^, 
which is negligible compared to the measured ve- 
locities. 

The radial velocity of Tycho-B was measured 
in the course of determining the stellar parame- 
ters of Tycho-B with the stellar parameter fitting 
package sfit (Jeffery et al. 2001). The sfit re- 
sult consistently gives Wheiio = —52 km s~^ for dif- 
ferent stellar parameters with an uncertainty of 
r-^ 2km s~^. 

In Table 3 we list all of the radial velocities 
both in a heliocentric frame and a local stan- 
dard of rest (LSR) frame. We will be referring 
to the heliocentric measurements henceforth. The 
listed uncertainty is the standard deviation of the 
radial-velocity measurement of all orders added in 
quadrature to the error of the radial- velocity stan- 
dards. 

In Figure 2 we compare the radial velocity of 
our sample stars to radial velocities of stars in 
the direction of Tycho's SNR using the Besangon 
Model (Robin et al. 2003). The distances as well as 
their imcertaintics arc taken from §3.6. The candi- 
dates' radial velocities are all typical for their dis- 
tances. Finally, we note that the measurement of 
Tycho-G is consistent with the results of WEK09 
and GH09. 

3.3. Rotational Velocity 

We have measured projected rotational veloci- 
ties {viot sin i) of all stars except Tycho-B in the 
fashion described by WEK09. We selected sev- 
eral unblended and strong (but not saturated) Fe I 
lines in the stellar spectra, and added them after 
shifting to the same wavelength and scaling to the 
same equivalent width (EW) . This was done to im- 
prove the S/N for the faint stars as well as provide 
consistency throughout all stars. 

As a reference we created three synthetic spec- 
tra for each star (one broadened only with the 
instrumental profile, the others with the instru- 
mental profile and a Uj-ot sini of 10 and 13km s^^, 
respectively) with the 2010 version of MOOG (Sne- 
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Table 2 

Proper Motions of Candidates 



Tycho 


a(J2000) 


5(J2000) 












l-M, 


V 


(Name) 


(hh;mm:ss.ss) 


(dd:mm:ss.ss) 


(mas yr -^) 


(mas yr ^) 


(mas yr ^) 


(") 


B 


0:25:19.97 


64:08:17.1 


-1.24 


0.56 


0.62 


0.64 


-1.17 


0.68 


4.86 


A 


0:25:19.73 


64:08:19.8 


-0.09 


-0.89 


1.17 


0.90 


-0.18 


-0.88 


6.21 


A2 


0:25:19.81 


64:08:20.0 


-U. ( i 


-O.DU 


U.oa 


U.D4 


-i.U < 


-o.oi 


6.58 


C 


0:25:20.38 


64:08:12.2 


-0.21 


-2.52 


0.65 


0.65 


-0.46 


-2.48 


6.66 


E 


0:25:18.28 


64:08:16.1 


2.04 


0.54 


0.66 


0.69 


2.09 


0.33 


7.60 


D 


0:25:20.62 


64:08:10.8 


-1.12 


-1.99 


1.01 


0.86 


-1.32 


-1.87 


8.60 


1 


0:25:16.66 


64:08:12.5 


-2.27 


-1.37 


1.60 


1.15 


-2.40 


-1.14 


18.00 


F 


0:25:17.09 


64:08:30.9 


-4.41 


0.20 


n 70 

U . ( u 


n 71 


4 "^7 




22.69 


J 


0:25:15.08 


64:08:05.9 


-2.40 


-0.25 


0.62 


0.62 


-2.42 


-0.00 


29.44 


G 


0:25:23.58 


64:08:01.9 


-2.50 


-4.22 


0.60 


0.60 


-z.yi 


-0.90 


29.87 


R 


0:25:15.51 


64:08:35.4 


0.28 


0.24 


0.89 


0.80 


0.30 


0.21 


33.23 


IN 


0:25:14. 7o 


c /( .no .00 1 

D4:Uo:zo.i 


1.18 


0.89 


0.86 


0.98 


1.27 


0.77 


00. DO 


u 


0:25:19.24 


64:07:37.9 


0.01 


-3.04 


0.73 


0.75 


-0.30 


-3.03 


36.06 


Q 


0:25:14.81 


64:08:34.2 


1.45 


3.07 


0.64 


0.72 


1.75 


2.91 


36.19 


T 


0:25:14.58 


64:07:55.0 


-3.85 


0.52 


0.72 


0.62 


-3.77 


0.91 


36.78 


K 


0:25:23.89 


64:08:39.3 


0.18 


0.17 


0.73 


0.69 


0.20 


0.15 


38.73 


L 


0:25:24.30 


64:08:40.5 


0.16 


-0.44 


0.75 


0.82 


0.11 


-0.45 


41.59 




u.^o. lo. i o 


64:08:34.4 


4.16 


0.58 


0.83 


0.84 


4.20 


0.15 


49 flQ 


9 




64:07:32.4 


74.85 


-4.43 


0.82 


0.83 


74.05 


-11.94 


4fi flQ 








Table 3 
















Radial Velocities of Candidates 








Tycho 


Date 


Vhelio 






a 










(Name) 


(dd/mm/yy) 


(dd:mm:ss.ss) 


(km s~ 


1) (k 


m s~^) 










Tycho-A 


09/09/06 


-36.79 


-28.5 





0.23 










Tycho-B 


09/09/06 


-52.70 


-44.41 


~ 2 










Tycho-C 


11/10/06 


-58.78 


-50.49 


0.75 










Tycho-D 


11/10/06 


-58.93 


-50.64 


0.78 










Tycho-E^ 


11/10/06 


-64.20 


-55.91 


0.27 










Tycho-G 


09/09/06 


-87.12 


-78.83 


0.25 










Tyclio-G 


ll/l()/()(i 


-(Sr.oi 


-79.2 


2 


0.78 











''There seems to be a discrepancy between RP04 and this work (measurement by RP04 t^LSR —26 km s ^), 
which might be a possible hint of a binary system. 
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/x« [mas/yr] 

Fig. 1. — The contours show the distribution of proper motions (68% and 95% probabihty) for all stars 
measured toward the Tycho SNR — excluding the named stars. We show the location of the candidate 
stars and their uncertainties on top of this distribution. Tycho-2 (called HP-1 in WEK09) is not shown in 
this figure as it is an extreme outlier with /i^ — 75 mas yr"-'^ and = —4.4 mas yr"-*^; it is also at a large 
distance from the remnant's geometric center (46"). We find, using the Besangon model as a proxy, that the 
contamination of this sample with foreground objects (less than 2kpc) is less than 15%. 



den 1973), using our derived temperature, surface 

gravity, and mctallicity. As input data to MOOG 
we used the Castelli & Kurucz (2004) atmospheric 
models and a line list from Kurucz & Bell (1995). 
We then applied the same process of line selection 
and adding as for the lines in the observed spectra. 

Figure 3 shows the comparison between the 
synthetic spectra with different rotational veloc- 
ities and the observed spectra. This comparison 
indicates that the stellar broadening (rotational, 
macroturbulcncc, etc.) is less than broadening due 
to the instrumental profile of 6 km s~^ for each 



star. We adopt 6 km s as an upper limit to the 

rotation for all stars. 

Due to its high temperature and rotation, we fit 
the rotational velocity for Tycho-B with the pro- 
gram SFIT (Jcffcry ct al. 2001, described in §3.4) as 
part of the overall fit for this star's stellar param- 
eters. We find Urotsinz = 171133 km s~^. While 
Tycho-B's rotation is very high compared to the 
other candidate stars, for stars of this temperature 
and surface gravity a high rotation is not unusual. 
In summary, other than Tycho-B, none of the stars 
show rotation which is measurable at this resolu- 




Distance [kpc] 



Fig. 2. — The fcontours indicate la, 2a, and 3a levels of the distance and radial velocity using the Besangon 
Model (Robin et al. 2003) with ~ 60,000 stars in the direction of SN 1572 (only including stars with 
10 < y < 20 mag and with a metallicity of [Fe/H] > -1 for the filled contours and [Fe/H] > -0.2 for 
the dashed contours). We have overplotted our candidate stars with error bars. One should note that the 
uncertainties in distance are a marginalized approximate of the error; the proper error surfaces can be seen 
in Figure 10. The vertical gray shade shows the error range for the distance to SNR 1572. 



tion. 

3.4. Stellar Parameters and Chemical 
Abundances 

The stellar parameters are presented in Table 5 
and were determined using a traditional spectro- 
scopic approach. Due to its high temperature, we 
measure the stellar parameters for Tycho-B by di- 
rect comparison to models in a separate procedure 
described later in this subsection. 

The first step in the spectroscopic analysis was 
to rectify the continuum. For each order, we fit 



the continuum, by eye, using a low-order poly- 
nomial function within the CONTINUUM task in 
IRAF. To help identify continuum regions in the 
program stars, wc made use of the Arcturus and 
solar spectra (Hinklc ct al. 2000). Consideration 
of the moderate S/N was a concern. For example, 
at these values of the S/N, we were mindful of not 
fitting the continuum to the highest points since 
it is likely that these values are noise rather than 
true continuum regions. 

Next, equivalent widths (EWs) for a set of 
Fe and Ni lines were measured using routines in 



8 




Relative Wavelength [km s ^] Relative Velocity [km s ^] 



Fig. 3. — The figures show the combination of iron line profiles after normalization to the same EW and 
compare them to synthetic line profiles created by MOOG. We convolved the synthetic lines first with a 
rotational kernel with three different values for rotation and then with the instrumental profile. All stars 
show rotation less than 6 km s~^, which is equal to the instrumental profile at this resolution. 



IRAF. The log gf values for the Fe I lines were 
from the laboratory measurements by the Oxford 
group (e.g., BlackweU et al. 1979a,b, 1980, 1986, 
1995, and references therein) and the Fe II lines 
were from the measurements by Biemont et al. 
(1991). For Ni, the log gf values were taken from 
the compilation by Reddy et al. (2003, henceforth 
ReddyOS) and Ramirez & Cohen (2002, henceforth 
RC02). While these EW measurement routines 
employ Gaussian fits in a semi-automated manner, 
we emphasize that all EWs were visually checked 
on at least two occasions. We also required that 
lines have an EW of at least 10 mA to avoid mea- 
suring noise and less than ^ 150 mA to avoid sat- 



urated lines with non-Gaussian profiles that may 
lie on the fiat part of the curve-of-growth. Table 8 
shows the EWs measured for the program stars. 
Missing values indicate that the line was not de- 
tected or that no reliable measurement could be 
obtained. In the following subsection, we consider 
in more detail the uncertainties that arise from 
continuum placement and EW errors. 

We used the 2011 version (Sobeck et al. 2011) of 
the local thermodynamic equilibrium (LTE) stel- 
lar line analysis program MOOG (Sneden 1973) 
and LTE model atmospheres from the Castelli & 
Kurucz (2003) grid to derive an abundance for a 
given line. The effective temperature, T^s, was 



9 



adjusted until the abundances from Fe I lines dis- 
played no trend as a function of lower excitation 
potential, %• The surface gravity, log g, was ad- 
justed until the abundances from Fc I and Fc II 
lines were in agreement. The microturbulent ve- 
locity, ^t, was adjusted until there was no trend 
between the abundances from the Fe I lines and 
the reduced EW, log (EW/A). This process was it- 
erated until self-consistent stellar parameters were 
obtained for each star. 

In our analysis, we explored stellar parame- 
ters at discrete values. For effective tempera- 
ture, we considered values at every 25 K (e.g., 
6000, 6025 K, etc.); for surface gravity, we consid- 
ered values at every 0.05 dex (e.g., 4.00, 4.05 dex, 
etc.); and for ^i, we considered values at every 
0.01km s-i (e.g., 1.70, 1.71km s-\ etc.). We 
assumed that excitation equilibrium was satisfied 
when the slope between log e(Fe I) and lower 
excitation potential (x) was < 0.004. We as- 
sumed that ionization equilibrium was achieved 
when log e(Fe I) - log e(Fe II) < 0.02 dex. The 
microturbulent velocity was set when the slope be- 
tween log e(Fe I) and reduced EW (log EW/A) was 
< 0.004. We found a unique solution for all pro- 
gram stars. We estimate that the internal uncer- 
tainties are typically Teff ± 100 K, \ogg ± 0.3 dex, 
and ± 0.3 km s~^. For further details regard- 
ing the derivation of stellar parameters, see Yong 
et al. (2008). 

The final iron measurements are the average of 
Fc I and Fe II, weighted by the number of lines 
measured for each species. We adopted the solar 
abundances of Asplund et al. (2009). 

In addition, we measured element abundance 
ratios for Ni via EW analysis and Li (only for 
Tycho-G) via spectrum synthesis (see Figure 6). 
For the Li spectrum synthesis, we used the Rcddy 
et al. (2002) line list in combination with MOOG 
and the Castelli & Kurucz (2003) model atmo- 
spheres. A non-LTE (NLTE) analysis (Lind et al. 
2009) of the Li abundances (A(Li)NLTE = 2.45) 
yields nearly the same result as the LTE abun- 
dance (A(Li) = 2.46). Abundances arc presented 
in Table 6. Tycho-B's abundances are not pre- 
sented in the table as they were measured in a 
different fashion. 

In summary, the inferred metallicities for all 
candidates show that the candidates are of roughly 
solar metallicity with the exception of the metal- 



poor Tycho-C. The range of metallicities spanned 
by the program stars is compatible with member- 
ship in the thin disk. Based on metallicity alone, 
we do not regard any of the program stars to be 
unusually metal-poor or metal-rich. Additionally, 
we have compared the [Ni/Fe] abundance ratio to 
a well-calibrated set of F- and G-dwarf abundances 
(Bensby et al. 2005), which we calibrated to the 
solar abundances of Asplund et al. (2009). Fig- 
ure 4 shows that all program stars are consistent 
with stars of similar metallicity. We do note that 
Tycho-C is a marginal outlier (perhaps Icr) with 
a low [Ni/Fc] abundance ratio, but do not regard 
this to be significant. To avoid selection effects 
we compared Tycho-C to a sample of giant stars 
(Luck & Heitcr 2007), which gives a similar result 
as the comparison with Bensby et al. (2005). 

Because Tycho-B has a temperature greater 
than 9000 K and is rapidly rotating, the process 
described above cannot be used to measure stel- 
lar parameters. Instead we used the program 
SFIT to match the HIRES spectrum to a grid 
of model spectra. To determine the stellar pa- 
rameters for Tycho-B we have used a model grid 
with [Fe/H] = -1.0, 8000 < Teff < 16,000 K. and 
7 < log g < 2. This low metallicity is suggested 
by the very weak Ca II K line and Mg II lines, but 
it is hard to measure. We cannot measure helium 
directly in this spectrum and thus adopt A'' (He) = 
0.1, as this is empirically a very common helium 
abundance in stars. 

This analysis resulted in Teff = 10, 000^2oo K, 
log g = 3.67 with slopeSlog g/dT^g = 0.27/500K"\ 
and rotational velocity t^rot sini = 171 km s~^ with 
slope dv.ot sin i/OT^a = -41/500 km s^^K^^ 
From qualitative analysis this object seems metal 
poor (e.g., in comparison to stars of similar stel- 
lar parameters hut solar metallicity), but its high 
rotation and temperature make it hard to deter- 
mine this parameter precisely. For the present, we 
assume [Fe/H] = —1.0 unless otherwise noted. 

In addition, using the high-resolution spectrum, 
we measured the EWs of several lines predicted 
to be strong in the Vienna Atomic Line Database 
(VALD; Kupka et al. 2000). The abundances were 
deduced from the EWs using a model atmosphere 
having T^s ^ 10,000 K, log g = 3.67, and [Fe/H] 
= -1.0 (see Table 4). 

One caveat regarding these abundances is the 
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Tycho-G 




[Fe/H] 

Fig. 4. The background gray error bars arc F- and G-dwarf abundances from Bensby et al. (2005). All 
candidate stars are consistent with that distribution. Tycho-C can be seen as an outlier, but it is a K-giant 
and its class is not represented in the underlying F- and G-dwarf distribution. 



use of EWs from single lines with large rota- 
tional broadening, since the effect of blending with 
nearby weak lines cannot be taken into account. A 
second is that these abundances invariably rely on 
the strongest lines, which are precisely those most 
susceptible to departures from LTE. Nevertheless, 
they do confirm the earlier impression that the star 
is metal-poor, and justify the adoption of [Fe/H] 
= -1.0 ±0.4. 

As a second approach to determine the stellar 

parameters of Tycho-B we used the low-resolution 
spectrum obtained with LRIS. The observation 
range of LRIS was chosen to be centered around 
the Balmcr jump, as this feature is sensitive to 
the surface gravity (Bessell 2007). We fitted the 



spectrum to a grid of model spectra (Munari et al. 
2005) using a spectrum-fitting tool described be- 
low. The final grid we used covered log g from 
3.5 to 4.5 in steps of 0.5 and effective tempera- 
ture from 9000 to 12,000 K in steps of 500 K. In 
addition, we expanded the grid by reddening the 
spectra with the pysynphot^ package. We also 
added diffuse interstellar bands (Bcals & Blanchet 
1937; Herbig 1966, 1967, 1975, 1995; Hibbins et al. 
1994; Jenniskens & Desert 1994; Wilson 1958) to 
the synthetic spectra, scaled with reddening. The 
included E{B — V) ranged from 0.5 to 1.3 mag in 



2 The PYSYNPHOT package is a product of the Space Tele- 
scope Science Institute, which is operated by AURA for 
NASA. 
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Table 4 
Tycho-B abundances 



Ion 

designation 


A 

(A) 


(A) 


e 

(dex) 


[X/H] 
(dex) 


de 
a log g 


de 

(K-1) 


Mg II 


4481.13+4481.33 


220 ± 15 


6.18+ .08 


-1.40 


0.08 


8 X 10-5 


Si II 


6347.1 


140 ±5 


6.96+ .18 


-0.59 


-0.02 


1 X 10"^ 


01 


7771.9+7774.2+7775.4 


460 ± 30 


8.43+ .10 


-0.58 


0.24 


-4 X 10-5 



steps of 0.2. We assumed a rotation of 171 km s 
in the grid (see §3.3). 

We used as a figure of merit in our fitting 
procedure. To find the best fit for Tycho-B we 
used the migrad algorithm provided by minuit 
and hnearly interpolated between the grid points 
using LinearNDInterpolator provided by the 
SciPY package. The fit of Tycho-B results in T^s 
= 10,570 K, log .9 = 4.05, [Fe/H] = -1.1, and 
E{B - V) = 0.85 mag. The model fits the syn- 
thetic spectrum poorly in the wavelength region 
3800-4280 A (see Figure 5). The adopted mixing- 
length parameter in one-dimensional (ID) model 
atmospheres, used to construct the spectral grid, 
influenc;c;s the fluxes in that region and affects the 
hydrogen line profiles. Heiter et al. (2002) and 
others show that a mixing length of 0.5, rather 
than 1.25 as used in the Kurucz/Munari grid, bet- 
ter fits the violet fiuxes and the H line profiles. 
Spectra using a mixing-length parameter of 0.5 are 
brighter in the ultraviolet, and the H(5, H7, and 
H/3 profiles give the same effective temperature 
as the Ha profiles. We have chosen, however, to 
fit the spectrum and ignore the problematic spec- 
tral region (3800-4280 A) to avoid a systematic 
error. This yields %s = 10, 722 K, log g = 4.13, 
[Fe/H] = -1.1, and E{B - V) = 0.86 mag. The 
differences are indicative of the size of system- 
atic errors in the model fits. We adopt the fit 
excluding the problematic wavelength region in 
the subsequent analysis. Exploring the complex 
search space, we estimate the uncertainties to be 
AT^ff = 200 K, A log g = 0.3, and A[Fe/H] = 0.5, 
and we note that the parameters are correlated. 
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1.0 




Tycho-B 
Best Fit 1 
Best Fit 2 



3500 



4000 4500 5000 
Wavelength A 



5500 



Fig. 5. — The plot shows the normahzed spectrum of Tycho-B with the fit which excluded the spectral region 
3800-4500 A (Best Fit 1) and the fit with the problematic region (Best Fit 2). The region is marked with a 
grey shade. 
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Table 5 
Stellar Parameters 



Tycho 




logs 


it 


[Fe/H] 




(K) 


(dex) 


(km s-^) 


(dex) 


A 


4975 


2.9 


1.20 


0.04 


B 


10722 


4.13 


2 


-1.1 


C 


4950 


2.9 


2.14 


-0.55 


E 


5825 


3.4 


1.82 


-0.13 


G 


6025 


4.0 


1.24 


-0.13 
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Table 6 

Chemical Abundances and Error Estimates for the Program Stars. 



Species 


l0gc(Jf)Q 


loge(X) 


[X/H] 


[X/Fc] 


a 












A [X/H] 


A[X/Fc] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) (8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


Star A 


Fe I 


7.50 


7.54 


0.04 




0.16 


40 0.03 


-0.08 


-0.01 


0.14 


0.02 


0.16 




Fe II 


7.50 


7.54 


0.04 




0.15 


14 0.04 


0.07 


-0.16 


0.12 


0.06 


0.22 




Ni I 


6.22 


6.33 


0.10 


0.06 


0.08 


37 0.01 


-0.05 


-0.04 


0.11 


0.03 


0.13 


0.04 


Star C 


Fe I 


7.50 


6.95 


-0.55 




0.23 


36 0.04 


-0.10 


-0.02 


0.07 


0.01 


0.13 




Fe II 


7.50 


6.96 


-0.54 




0.27 


5 0.12 


0.05 


-0.16 


0.05 


0.06 


0.22 




Ni I 


6.22 


5.52 


-0.70 


-0.15 


0.21 


20 0.05 


-0.06 


-0.04 


0.07 


0.03 


0.10 


0.05 


Star E 


Fe I 


7.50 


7.37 


-0.13 




0.22 


40 0.03 


-0.09 


0.02 


0.11 


0.00 


0.14 




Fe II 


7.50 


7.38 


-0.12 




0.16 


7 0.06 


0.03 


-0.14 


0.08 


0.03 


0.18 




Ni I 


6.22 


6.1 


-0.12 


0.01 


0.17 


8 0.06 


-0.07 


0.00 


0.09 


0.00 


0.11 


0.04 


Star G 


Fe I 


7.50 


7.37 


-0.13 




0.18 


69 0.02 


-0.09 


0.03 


0.09 


0.01 


0.13 




Fe II 


7.50 


7.35 


-0.15 




0.18 


16 0.05 


0.01 


-0.12 


0.07 


0.03 


0.15 




Ni I 


6.22 


6.16 


-0.06 


0.07 


0.14 


18 0.03 


-0.07 


0.01 


0.05 


0.01 


0.09 


0.04 


Li 




2.46 








1 0.05* 


-0.07 


0.01 


0.00 


0.02 


0.09 





^This represents the uncertainty in the fitting of the Une. 



15 



3.5. Tycho-G: A Detailed Comparison 
with GH09 

GH09 suggested that Tycho-G is a plausible 
donor star, with the primary evidence consisting of 
an unusually high Ni abundance and a high space 
velocity (radial velocity and proper motion). In 
this subsection, we focus on the Ni abundance, 
and we refer the reader to §3.1, 3.2, and 4 on the 
proper motion and radial velocity. 

The measured values are [Ni/Fe] = 0.16 ± 0.04 
and 0.07 ± 0.04 for GH09 and this study, respec- 
tively, from the same HIRES spectra. The mag- 
nitude of the difference is 0.09 dex, and it is sig- 
nificant at the ~ 1.5(7 level. While our [Ni/Fe] 
ratio in Tycho-G is lower than that measured by 
GH09, our value does not represent a substantial 
revision given the measurement uncertainties in- 
volved. Nevertheless, our [Ni/Fe] measurement 
and comparison with the literature do not support 
an unusually high Ni abundance, and we conclude 
that Tycho-G does not show any obvious chem- 
ical signature that one may seek to attribute to 
a supernova companion star. In order to identify 
the origin of the difference in [Ni/Fe] ratios, we 
now compare our stellar parameters and chemical 
abundances to those of GH09. 

Both studies determined stellar parameters and 
chemical abundances in a similar manner, from a 
standard spectroscopic EW analysis using ID LTE 
Kurucz model atmospheres and the MOOG stel- 
lar line analysis software. Our analysis employed 
more recent versions of both tools. The first test 
we can perform is to use the GH09 line list and 
stellar parameters but with our tools — namely, 
the 2011 version of MOOG (Sobeck et al. 2011; 
Sneden 1973) and the Castelh & Kurucz (2003) 
model atmospheres. Adopting this approach, we 
obtain log e(Fe I) = 7.38 (ct = 0.13), log e(Fe II) 
= 7.42 (cr = 0.10), and log e(Ni I) = 6.33 {a = 
0.19). These values are in very good agreement 
with those of GH09, who obtained log e(Fe I) = 
7.42 (cr = 0.12), log e(Fe II) = 7.42 (cr = 0.10), 
and log e(Ni I) = 6.36 {a = 0.19). Thus, we ar- 
gue that any abundance differences (for Fe and Ni) 
between the two studies, exceeding the ^0.04 dex 
level, cannot be attributed to differences in the 
model-atmosphere grid and/or line-analysis soft- 
ware. 

Our stellar parameters (reff= 5900 ± lOOK, 



log g = 3.85 ±0.30, [Fe/H] = -0.05 ±0.09) are in 
good agreement with those of GH09 (Teff= 6000 ± 
100, log g = 4.00 ± 0.30, [Fe/H] = -0.13 ± 0.13). 
The second test we can perform is to determine 
chemical abundances using (i) the GH09 stellar 
parameters but with our line list and (ii) our stel- 
lar parameters and line list. On comparing case 
(ii) minus case (i), we find A log e(Fe I) = 0.10, 
A log e(Fe II) = 0.02, and A log e(Ni I) = 0.08. 
Adopting the same solar abundances and method 
for determining the average [Fe/H] value (average 
of Fe I and Fe II weighted by the number of lines) 
as in the present study, we find A [Ni/Fe] = 0.00. 
We argue that while there are abundance differ- 
ences for log e{X) at the ~ 0.10 dex level, the 
[Ni/Fe] ratio remains unchanged, and therefore 
any differences in the [Ni/Fe] ratio between the 
two studies cannot be attributed to differences in 
the adopted stellar parameters. 

The solar abundances for Fe and Ni differ be- 
tween the two studies. GH09 adopt 7.47 and 6.25 
for Fe and Ni, respectively, while we use 7.50 and 
6.22 (from Asplund et al. 2009). Had we used the 
GH09 solar abundances, we would have obtained 
a ratio [Ni/Fe] = 0.01. Therefore, the different so- 
lar abundances adopted by the two studies only 
serve to decrease the discrepancy in the [Ni/Fe] 
ratio — that is, any difference in [Ni/Fe] cannot 
be attributed to the solar abundances. 

The next series of comparisons we can perform 
concern the line lists. We measured Fe and Ni 
abundances using the GH09 line list but with our 
stellar parameters and find log e(Fe I) = 7.42 (a = 
0.12), log e(Fe II) = 7.42 (cr = 0.10), and log e(Ni 
I) = 6.36 (cr = 0.19). Table 7 gives a comparison 
of all tests performed. 

Adopting the same approach as before, regard- 
ing the solar abundances and metallicity, yields a 
ratio [Ni/Fe] = 0.22, a value that exceeds both 
our measurement and that of GH09. We therefore 
speculate that the difference in [Ni/Fe] between 
the two studies is driven primarily by differences 
in the line list. In particular, we note that while 
the Fe I and Fe II abundances are in fair agree- 
ment with our value and GH09, it is the Ni abun- 
dance, log e(Ni), that shows a large difference be- 
tween the two studies: 6.16 ± 0.09 and 6.33 ± 0.10 
for this study and GH09, respectively. Although 
the magnitude of this difference may appear large, 
0.17 dex, it is significant only at the 1.3cr level. 
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On comparing the line lists between the two 
studies, we find 3, 2, and 8 lines in common for 
Fe I, Fe II, and Ni, respectively. For these three 
species, the log gf values are on the same scale 
with differences (this study minus GH09) of —0.04 
(c7 = 0.07), -0.03 (fT = 0.04), and -0.01 [a = 0.03) 
for Fe I, Fe II, and Ni, respectively. Although the 
comparison sample is small, there is no clear evi- 
dence for any large systematic difference in log gf 
values that could explain the differing log e(Ni) or 
[Ni/Fe] values. 

For the lines in common, our EWs are, on av- 
erage, lower than those of GHG9 by 5.7 mA (cr = 
8.0mA), 5.6mA {a = 5.4mA), and 12.7mA {a 
= 6.9 mA) for Fe I, Fe II, and Ni, respectively. 
The most intriguing aspect of this comparison is 
that the Ni lines show the greatest discrepancy. In 
light of the EW differences for Fe I and Fe II, we 
may naively have expected the Ni EWs to show 
an offset of ~6mA rather than a 12.7mA offset. 
Indeed, differences in the Ni EWs appear to be the 
primary reason for the difference in the derived Ni 
abundances between the two studies. 

In Figure 7, we plot our EWs and the GH09 
EWs, for the 8 Ni lines in common. To estimate 
the uncertainties in our EWs, we use the Cayrel 
(1988) formula which considers the measurement 
uncertainty due to the line strength, S/N, and 
spectral resolution. Uncertainty in the continuum 
placement is not included in the Cayrel (1988) for- 
mula. 

As noted in the previous subsection, we regard 
continuum placement to be an additional source of 
uncertainty in the EW measurements. To quan- 
tify this uncertainty, we use the DAOSPEC pro- 
gram which fits the continuum and measures EWs 
(Stetson & Pancino 2008). Using DAOSPEC, we 
remeasure the Ni EWs using four different contin- 
uum fitting criteria: (i) adopting our continuum 
placement, and using a (ii) third-order, (iii) fifth- 
order, and (iv) ninth-order polynomial to refit our 
continuum-rectified spectra. For a given line, we 
compute the dispersion in the EW measurements 
from the four different methods for continuum fit- 
ting and adopt this value as being representative 
of the EW uncertainties due to continuum recti- 
fication. We then add this value, in quadrature, 
to the uncertainty using the Cayrel (1988) value, 
noting that the latter value dominates the total 
EW error budget (see Table 9). 



To establish whether these EW imcertainties 
are valid, we first identify the set of Ni EWs which 
produce our mean [Ni/Fe] ratio. That is, every 
line in this set of "ideal" EWs produces log e(Ni) 
= 6.16, i.e., [Ni/Fe] = 0.07. We then added to 
each of these ideal EWs a random number drawn 
from a normal distribution of width correspond- 
ing to our estimate of the EW uncertainty. We 
repeated this process for each Ni line, computed 
Ni abundances for this new set of lines, and mea- 
sured the abundance dispersion. We repeated this 
process for 1,000 new random samples. The av- 
erage dispersion in Ni abundance is 0.17 dex (a 
= 0.06 dex), and this average value agrees well 
with our observed dispersion of 0.14 dex. There- 
fore, we are confident that our EW measurement 
uncertainties are realistic, since this Monte Carlo 
analysis verifies that these uncertainties reproduce 
our observed abundance dispersion. 

An additional test is to measure EWs from our 
spectra for all Fe and Ni lines measured by GH09. 
As with our EWs, all lines were manually checked. 
For Fe I, we measured 27 lines and found a mean 
difference (this study minus GH09) of-1.9mA ± 
1.2 (cr = 6.0). For Fe II, we measured 8 lines and 
found a mean difference of —4.6 mA ± 2.8 (cr = 
7.8). For Ni, we measured 18 lines and found a 
mean difference of -8.7mA ± 2.0 (cr = 8.4). This 
comparison confirms that our EWs are systemat- 
ically lower than those of GH09 and that the Ni 
lines, in particular, show the largest discrepancy. 
Indeed, the average difference in Ni EWs is 4 times 
larger than the average difference in Fe I EWs. 
While continuum normalization could potentially 
explain these differences, these Ni lines lie in spec- 
tral regions similar to those of the Fe lines, so we 
would expect the differences in EWs for Fe and Ni 
to behave similarly. 

We note in our line selection that we reject 5, 
2, and 4 lines of Fe I, Fe II, and Ni (respectively) 
that were measured by GH09. These lines were in 
our opinion blended and/or in regions where the 
local continuum was poorly defined. 

We return now to the eight Ni lines in common, 
noting that (i) for 7 of the 8 lines, our EWs are 
smaller than those of GH09, (ii) for 7 of the 8 
lines, the difference in EWs exceeds Icr, and for 
all 7 lines, the difference shows the same "sign," 
and (iii) for 4 of the 8 lines, the difference in EWs 
exceeds 2a, and for all 4 lines, the difference shows 
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the same "sign." 

Finally, for the eight Ni lines in common with 
GH09, we plot our normalized spectra along with 
spectrum syntheses (see Figures 8 and 9). The 
main points to take from these figures are the lo- 
cation of the continuum and how well the spec- 
trum syntheses fit the lincis for the abundances we 
measure. We note that our abundances were de- 
termined from EW analysis rather than spectrum 
synthesis. Nevertheless, had we relied solclj^ upon 
spectrum synthesis, we would have obtained es- 
sentially identical results. A systematic increase 
in log e(Ni) of 0.17 dcx or in [Ni/Fe] of 0.09 dex, 
as measured by GH09, is not supported by these 
spectrum syntheses. 

The main conclusions we draw from this com- 
parison arc (i) abundance differences between the 
two studies cannot be attributed to the differ- 
ent versions of model atmospheres and spectrum 
synthesis software; (ii) the [Ni/Fc] ratio remains 
unchanged when using our line list but with ei- 
ther the GH09 stellar parameters or our stellar 
parameters; (iii) differences in [Ni/Fe] cannot be 
attributed to the adopted solar abundances; (iv) 
although the set of lines in common between the 
two analyses is small, there are no large systematic 
differences in the log gf values that could explain 
the discrepancy in Ni abundances; (v) for Fe I and 
Fe II, our EWs are systematically lower than those 
of GH09 by - 6 mA, and our Ni EWs are system- 
atically lower by ^ 12 mA; and (vi) our EW un- 
certainties for Ni are consistent with the observed 
dispersion in Ni abundance. 

As noted above, while our measured [Ni/Fe] 
value does not represent a substantial revision of 
the GH09 value, our Ni abundance is not unusual 
with respect to field stars at the same metallic- 
ity. Nevertheless, we welcome further analyses of 
this star, preferably conducted with higher-quality 
spectra. 



18 



Table 7 

Comparison of Ni measurement 



Species 

(1) 


log eiX) 
(2) 


[X/H]* a 
(3) 


This work 


Fe I 
Fe II 
Ni I 


7.37 
7.35 
6.16 


-0.13 0.18 
-0.15 0.18 
-0.06 0.14 


Measurements from GH09 


Fe I 
Fe II 
Ni I 


7.42 
7.42 
6.36 


-0.08 0.10 
-0.08 0.12 
0.16 0.19 


This work using GH09's line hst and stellar parameters 


Fe I 
Fe II 
Ni I 


7.38 
7.42 
6.33 


-0.12 0.13 
-0.08 0.10 
0.11 0.19 


This work 


using GH09's line list and this work's stellar parameters 


Fe I 
Fe II 
Ni I 


7.42 
7.42 
6.36 


-0.08 0.12 
-0.08 0.10 
0.14 0.19 



^Using solar values from Asplund et al. (2009). 
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Table 8 

Line List and Equivalent Width Measurements 



velength (A) 


Species 


X (eV) 


logs/ 


Star A 


Star C 


Star E 


Star G 


Source 












EW ( 


mA) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


4602.00 


26.0 


1.607 


— 3.154 


114.8 








Oxford 


4733.59 


26.0 


1.484 


— 2.987 






102.3 


85.4 


Oxford 


4802.88 


26.0 


3.692 


— 1.531 


91.0 


96.0 




64.8 


Oxford 


4848.88 


26.0 


2.277 


— 3.154 




85.9 






Oxford 


4930.31 


26.0 


3.957 


— 1.264 








57.8 


Oxford 


4962.57 


26.0 


4.175 


— 1.199 


77.5 








Oxford 


5014.94 


26.0 


3.940 


— 0.320 








106.8 


Oxford 


5044.21 


26.0 


2.849 


— 2.034 


109.1 






67.2 


Oxford 


5049.82 


26.0 


2.277 


— 1.372 








117.8 


Oxford 


5054.64 


26.0 


3.637 


— 1.938 


70.8 






30.4 


Oxford 


5083.34 


26.0 


0.957 


— 2.958 








96.2 


Oxford 


5141.74 


26.0 


2.422 


— 2.001 


123.9 






85.4 


Oxford 


5151.91 


26.0 


1.010 


— 3.322 








85.0 


Oxford 


5166.28 


26.0 


0.000 


—4.195 






111.6 




Oxford 


5194.94 


26.0 


1.556 


— 2.090 








113.8 


Oxford 


5198.71 


26.0 


2.221 


— 2.135 








92.1 


Oxford 


5217.39 


26.0 


3.209 


— 1.179 






123.7 


81.9 


Oxford 


5223.19 


26.0 


3.632 


— 1.800 


58.5 




33.9 


27.7 


Oxford 


5225.52 


26.0 


0.110 


—4.789 






82.3 


65.6 


Oxford 


5242.49 


26.0 


3.632 


— 0.984 


113.9 




110.6 


69.8 


Oxford 


5247.05 


26.0 


0.087 


—4.946 






76.6 


55.5 


Oxford 


5250.21 


26.0 


0.121 


—4.938 








59.9 


Oxford 


5253.46 


26.0 


3.281 


— 1.630 








62.4 


Oxford 


5412.80 


26.0 


4.431 


— 1.783 


38.2 








Oxford 


5491.84 


26.0 


4.183 


— 2.253 


36.9 








Oxford 


5497.52 


26.0 


1.010 


— 2.849 






129.3 


102.1 


Oxford 


5501.46 


26.0 


0.957 


— 3.063 








90.9 


Oxford 


5506.78 


26.0 


0.989 


— 2.797 








95.7 


Oxford 


5525.54 


26.0 


4.227 


— 1.149 








38.4 


Oxford 


5569.62 


26.0 


3.414 


— 0.544 






123.1 


113.4 


Oxford 


5586.76 


26.0 


3.366 


— 0.161 






142.8 




Oxford 


5600.23 


26.0 


4.257 


— 1.486 








44.2 


Oxford 


5618.63 


26.0 


4.206 


— 1.292 


80.1 






43.3 


Oxford 


5661.35 


26.0 


4.281 


— 1.822 


48.2 






18.0 


Oxford 


5662.51 


26.0 


4.175 


-0.590 


115.6 




79.1 


69.5 


Oxford 


5701.55 


26.0 


2.557 


-2.216 






104.7 




Oxford 


5705.47 


26.0 


4.298 


-1.421 


78.9 


69.5 






Oxford 


5741.85 


26.0 


4.253 


-1.689 


58.0 








Oxford 


5753.12 


26.0 


4.257 


-0.705 








52.0 


Oxford 


5775.08 


26.0 


4.217 


-1.314 


82.3 


73.0 




34.5 


Oxford 


5778.45 


26.0 


2.586 


-3.481 


55.5 


58.5 






Oxford 


5816.37 


26.0 


4.545 


-0.618 








66.5 


Oxford 


5855.09 


26.0 


4.604 


-1.547 


43.1 




17.1 




Oxford 


5909.97 


26.0 


3.209 


-2.643 




44.9 




36.5 


Oxford 


5916.25 


26.0 


2.452 


-2.994 




73.6 




44.1 


Oxford 


5956.69 


26.0 


0.858 


-4.608 


107.2 


101.5 


59.7 


34.9 


Oxford 


6012.21 


26.0 


2.221 


-4.073 


72.0 








Oxford 


6027.05 


26.0 


4.073 


-1.106 


100.3 


79.9 




63.3 


Oxford 


6065.48 


26.0 


2.607 


-1.530 






106.2 




Oxford 


6082.71 


26.0 


2.221 


-3.573 




61.6 






Oxford 


6120.24 


26.0 


0.914 


-5.970 


33.9 








Oxford 


6136.62 


26.0 


2.452 


-1.400 






145.6 




Oxford 


6136.99 


26.0 


2.196 


-2.950 






71.7 




Oxford 


6151.62 


26.0 


2.174 


-3.299 


91.4 


81.8 


33.7 


36.6 


Oxford 


6165.36 


26.0 


4.140 


-1.490 


73.5 




54.0 




Oxford 


6173.34 


26.0 


2.221 


-2.880 


115.1 






52.8 


Oxford 
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Table 8 — Continued 



Wavelength (A) 


Species 


X (cV) 


logg/ 


Star A Star C Star E Star G 


Source 


(1) 


(2) 


(3) 


(4) 


EW (mA) 
(5) (6) (7) (8) 


(9) 



6180.20 


26.0 


2.725 


— 2.637 




61.0 




46.8 


Oxford 


6200.31 


26.0 


2.607 


— 2.437 




109.6 




66.6 


Oxford 


6219.28 


26.0 


2.196 


— 2.433 




134.4 




75.2 


Oxford 


6229.23 


26.0 


2.843 


— 2.846 


82.3 




61.5 




Oxford 


6230.73 


26.0 


2.557 


— 1.281 








109.9 


Oxford 


6232.64 


26.0 


3.651 


— 1.283 


114.5 


129.1 


106.6 


65.3 


Oxford 


6246.32 


26.0 


3.600 


— 0.894 








106.2 


Oxlord 


6252.55 


26.0 


2.402 


— 1.687 








100.8 


Oxlord 


6265.13 


26.0 


2.174 


— 2.550 








75.0 


Oxford 


6270.22 


26.0 


2.856 


— 2.505 




77.1 


70.2 




Oxford 


6297.79 


26.0 


2.221 


— 2.740 


120.6 


115.1 


93.5 


61.7 


Oxford 


6301.50 


26.0 


3.651 


— 0.766 






114.0 


110.9 


Oxford 


6322.69 


26.0 


2.586 


— 2.426 


120.9 


119.1 






Oxford 


6335.33 


26.0 


2.196 


— 2.194 






116.4 


75.0 


Oxford 


6336.82 


26.0 


3.684 


— 0.916 








69.1 


Oxford 


6344.15 


26.0 


2.431 


— 2.923 








48.4 


Oxford 


6355.03 


26.0 


2.843 


— 2.403 




101.7 


93.4 


64.4 


Oxford 


6393.60 


26.0 


2.431 


— 1.469 






129.4 


105.5 


Oxford 


6408.02 


26.0 


3.684 


— 1.066 




127.1 


87.4 


80.7 


Oxford 


6411.65 


26.0 


3.651 


— 0.734 








118.0 


Oxford 


6430.84 


26.0 


2.174 


— 2.006 






135.0 


102.5 


Oxford 


6481.87 


26.0 


2.277 


— 2.984 


113.2 






50.3 


Oxford 


6498.94 


26.0 


0.957 


— 4.687 




121.1 




35.3 


Oxford 


6574.23 


26.0 


0.989 


— 5.004 


84.3 


88.5 




25.3 


Oxford 


6575.02 


26.0 


2.586 


— 2.727 


108.2 


115.3 




51.0 


Oxford 


6592.91 


26.0 


2.725 


— 1.490 






121.7 


104.1 


Oxford 


6593.87 


26.0 


2.431 


— 2.422 






99.1 


75.4 


Oxford 


6609.11 


26.0 


2.557 


— 2.692 




104.9 


54.7 


53.1 


Oxford 


6648.08 


26.0 


1.010 


— 5.918 


48.2 








Oxford 


6677.99 


26.0 


2.690 


— 1.435 






142.4 


98.3 


Oxford 


6699.16 


26.0 


4.590 


— 2.170 


24.6 


15.8 






Oxford 


6739.52 


26.0 


1.556 


— 4.823 


48.4 


53.3 






Oxford 


6750.15 


26.0 


2.422 


-2.621 


120.0 


96.6 


70.2 


60.0 


Oxford 


6752.70 


26.0 


4.635 


-1.273 








33.8 


Oxford 


6810.26 


26.0 


4.603 


-1.003 


73.7 






37.5 


Oxford 


6837.02 


26.0 


4.590 


-1.756 




19.0 


21.8 




Oxford 


7112.17 


26.0 


2.988 


-3.044 


86.4 


48.9 


40.7 




Oxford 


7223.66 


26.0 


3.015 


-2.269 




87.5 


49.9 


35.0 


Oxford 


7401.69 


26.0 


4.183 


-1.664 




38.8 


30.4 




Oxford 


7710.36 


26.0 


4.217 


-1.129 




66.5 


56.2 




Oxford 


7723.20 


26.0 


2.277 


-3.617 




86.2 


37.3 




Oxford 


7912.86 


26.0 


0.858 


-4.848 




97.1 




21.7 


Oxford 


7941.09 


26.0 


3.271 


-2.331 


76.6 


42.5 




23.9 


Oxford 


8075.15 


26.0 


0.914 


-5.088 


105.4 


118.7 






Oxford 


4491.40 


26.1 


2.853 


-2.684 






106.6 




Biemont93 


4508.29 


26.1 


2.853 


-2.312 


104.1 


106.2 




108.0 


Biemont93 


4620.52 


26.1 


2.826 


-3.079 


75.7 






69.7 


Biemont93 


4635.32 


26.1 


5.952 


-1.275 






48.1 


23.9 


Biemont93 


4993.36 


26.1 


2.805 


-3.485 


58.4 








Biemont93 


5100.66 


26.1 


2.805 


-4.135 








42.1 


Biemont93 


5132.67 


26.1 


2.805 


-3.901 


41.5 








BiemontQS 


5197.58 


26.1 


3.228 


-2.233 


101.2 






86.8 


Bicmont93 


5234.62 


26.1 


3.219 


-2.151 


105.4 






93.4 


Bicmont93 


5414.07 


26.1 


3.219 


-3.750 


45.8 








Biemont93 


5425.26 


26.1 


3.197 


-3.372 


58.1 


63.6 






Biemont93 


5991.38 


26.1 


3.150 


-3.557 




30.0 


54.4 


23.5 


Biemont93 


6084.11 


26.1 


3.197 


-3.808 


37.3 








Biemont93 
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Table 8 — Continued 



Wavelength (A) 


Species 


X (oV) 


log gf 


Star A 


Star C 


Star E 


Star G 


Source 












EW ( 


T\ 

mA) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


u j-iy . 


26 1 


3 886 


— 2 724 


48 1 




53 2 


38 2 


DltJlIKJlltyO 


6239 95 


26 1 


3 886 


3 439 








25 8 


Biern.orit93 


6247 56 


26 1 


3 889 


— 2 329 


60 4 




74 7 


68 7 


BicrQoiit93 


OOO c) ."^D 


26 1 


2 889 


— 4 253 


38 5 


28 2 


41 7 


24 1 


Bicrn.orit93 


DOOO . 1 Z 


26 1 


5 548 


2 271 








12 4 


Biern.oii.t93 




26 1 


2 889 


2 708 


62 8 


39 






Biem.ont93 


O'tOD.OO 


26 1 


3 900 


— 2 075 








75 7 


Bi6monti93 


fif^l Oft 
DOlD.Uo 


26 1 


9 S3SQ 

z.ooy 


— O.40U 










Biemont93 


7222 39 


26 1 


3 886 


3 295 
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Bicm.orit93 


7711 72 


26 1 


3 900 


2 543 


68 6 




85 7 


58 9 


B i ciiio lit 9 3 


'^089 '^'^ 
OUoZ . oO 


28 


3 660 


Q 59Q 


87 3 




69 2 


65 5 


ii.CU.iJ.j'UO 


OUOO.O'i 


Zo.U 


O.oDU 


-j^ 040 


Oo.z 






■^n 9 
ou.z 


iXCQQyUO 


ouoo.yu 


Zo.U 


O.DoU 


-j^ 240 


0O.4I: 








ixCQQyUo 




Zo.U 


O.OOU 


1 n7n 

— 1 .U ( u 


48 3 








xxCQOyUo 


5115 40 


28 


3 830 


Q 280 


89 1 








xxcQciyuo 




28 


4 100 


— 470 


75 5 








RC02 


/ ftO.OO 


98 n 

^o.U 


1 .OoU 


— O.ZDU 


7/1 7 


62 2 










28 


3 940 


990 


17 1 








RC02 


0041: ( .Ui 


Zo.U 


i .DoU 


/t 1 n 

t5.4:lU 




oy .0 






IxCQQyUo 


6007 31 


28 


1 680 


3 340 


66 


37 9 






RC02 




28 


4 230 


]^ 070 


40 5 








RC02 


DUoD . zo 


28 


4 260 


Q 520 


Di.y 








iXV_/UZ 


D 1 Uo . IZ 


28 


1 680 


2 440 


1118 


94 6 




[^■^ fi 

OO.D 


IXV-'UZ 


Dlll.Uo 


Zo.U 




n SI n 

— U.olU 


fin Q 
DU.y 


22 5 






rteQayuo 


6130 14 


28 


4 270 


— 940 


43 9 








r\.eQQyuo 


6175 37 


28 


4 090 


— 550 


69 6 


55 7 




46 8 


rteciciyuo 


D -L f O . oZ 


OS n 

Zo.U 


4 090 


— U.ZDU 


Si^ s 


pi7 /I 


7 

OD. ( 


t^7 /I 


xxeQQyuo 


6177 25 


28 


1 830 


3 510 


50 3 








rvc; <-n_i y u o 


6186 71 


28 


4 100 


Q 970 


52 7 






21 6 


RC02 


U ^ 4: . U J- 


28 


4 090 


110 


46 2 








rv t; <-n_i y u o 


6322.17 


28.0 


4.150 


-1.170 


39.5 








RC02 


6370.35 


28.0 


3.540 


-1.940 


38.4 








RC02 


6378.26 


28.0 


4.150 


-0.830 


56.1 


40.9 






Reddy03 


6482.80 


28.0 


1.930 


-2.630 


88.8 


64.4 




24.1 


RC02 


6598.60 


28.0 


4.230 


-0.980 


44.8 


24.9 






RC02 


6635.12 


28.0 


4.420 


-0.830 


45.1 








RC02 


6643.64 


28.0 


1.680 


-2.030 


151.0 


127.3 


95.1 


71.8 


Reddy03 


6767.77 


28.0 


1.830 


-2.170 


127.1 


129.1 




66.9 


RC02 


6772.32 


28.0 


3.660 


-0.970 


78.0 








Reddy03 


6842.04 


28.0 


3.660 


-1.470 


58.1 






23.9 


RC02 


7030.01 


28.0 


3.540 


-1.730 


53.6 








RC02 


7122.20 


28.0 


3.540 


0.050 


146.4 


120.8 






RC02 


7261.92 


28.0 


1.950 


-2.700 








38.1 


RC02 


7327.65 


28.0 


3.800 


-1.770 


35.1 






13.3 


RC02 


7409.35 


28.0 


3.800 


-0.100 




83.6 






RC02 


7414.50 


28.0 


1.990 


-2.570 




83.6 


48.0 


40.4 


RC02 


7422.28 


28.0 


3.630 


-0.130 


139.5 


101.4 


87.4 




RC02 


7574.05 


28.0 


3.830 


-0.580 


91.2 


39.5 




41.8 


RC02 


7748.89 


28.0 


3.700 


-0.380 


117.7 


85.7 


95.7 


70.3 


Reddy03 


7788.93 


28.0 


1.950 


-2.420 




123.4 


82.8 




RC02 


7797.59 


28.0 


3.900 


-0.350 


106.4 


84.0 


81.1 


61.7 


Reddy03 


7917.44 


28.0 


3.740 


-1.500 








14.8 


RC02 



Note.— Oxford = Blackwell et al. (1979b, a, 1980, 1986, 1995) 
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Table 9 

Equivalent Width Uncertainties for Ni in Star G 



Wavelength (A) 


EW 








(1) 


(2) 


(3) 


(4) 


(5) 


5082.35 


55.1 


5.0 


1.2 


5.1 


5088.54 


25.7 


5.0 


0.8 


5.0 


6086.28 


33.5 


6.0 


1.3 


6.1 


6108.12 


58.1 


5.8 


1.1 


5.9 


6175.37 


39.8 


6.4 


1.1 


6.5 


6176.82 


54.4 


6.0 


1.0 


6.0 


6186.71 


21.1 


5.2 


1.9 


5.6 


6482.80 


37.6 


6.0 


0.7 


6.0 


6643.64 


79.9 


6.2 


1.1 


6.3 


6767.77 


66.9 


6.8 


2.4 


7.2 


6842.04 


18.7 


6.3 


1.7 


6.6 


7261.92 


35.6 


6.7 


1.3 


6.9 


7327.65 


8.4 


6.4 


1.5 


6.6 


7414.50 


40.9 


7.0 


0.2 


7.0 


7574.05 


53.8 


6.4 


0.8 


6.5 


7748.89 


71.3 


7.4 


0.7 


7.5 


7797.59 


63.3 


7.7 


1.4 


7.8 


7917.44 


16.8 


7.5 


0.4 


7.5 



^This is the error from the Cayrel (1988) for- 
mula. 

^This is the error due to continuum placement 
(see text for details). 

■^This is the total error oljiaincd by adding 
columns (4) and (5) in quadrature. 
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3.6. Distances 

To measure the distance to the candidate stars 
we used colors and absolute magnitudes from 
isochrones by Pietrinferni et al. (2004). We used 
the MIGRAD algorithm (James & Roos 1975) to 
find close matches of the measured values to 
TefT-log g isochrones by varying the age of the 
isochrone. Subsequently we calculated E{B — V) 
using the isochrone's color, and we extracted a 
mass from the isochrone. The results can be seen 
in Table 10. To estimate the uncertainties in all 
distances, reddenings, and masses, we employed 
the Monte-Carlo method with 10,000 samples of 
effective temperature, surface gravity, metallicity, 
B magnitude, and V magnitude (see Figure 10). 
Errors included in Table 10 are the standard de- 
viations of the Monte-Carlo sample. The data 
show that all stars are compatible with the dis- 
tance of the remnant. This is not unexpected, as 
the uncertainties of the measurements in stellar 
parameters are relatively large. 

4. Discussion 

In our sample of six stars we find no star that 
shows characteristics which strongly indicate that 
it might be the donor star of SN 1572. On the 
other hand, it is difficult to absolutely rule out any 
particular star, if one is able to invoke improbable 
post-explosion evolutionary scenarios. 

Tycho-A is a metal-rich giant, and it seems 
likely to be a foreground star. Its principal re- 
deeming feature as a donor-star candidate is that 
it is located in the geometric center of the rem- 
nant, and that it has a relatively low surface grav- 
ity. Tycho-A shows a very low spatial motion, 
which is consistent with a giant-donor-star sce- 
nario, although its lack of rotation is in conflict 
with a donor-star scenario. Taking all measure- 
ments into account, we regard Tycho-A to be a 
very weak candidate (although a wind accretion 
scenario might still work). 

Tycho-B's high temperature, position at the 
center of the remnant, high rotational velocity, 
and unusual chemical abundance make it the most 
unusual candidate in the remnant's center. De- 
spite the a posteriori unlikely discovery of such a 
star in the remnant's center, Tycho-B's high rota- 
tional velocity coupled with its low spatial velocity 
seem to be in conflict with any viable donor-star 



scenario. These scenarios predict that the donor 
star will tidally couple to the white dwarf before 
explosion, causing the rotation and spatial mo- 
tion to be correlated post explosion (as discussed 
in WEK09). The large rotation seen in Tycho-B 
should be accompanied by a large spatial motion, 
which is ruled out by the observations presented 
here, a problem we are unable to reconcile with 
Tycho-B being the donor star. However, Tycho- 
B docs show some unusual abimdances, which we 
will scrutinize in future studies. 

Tycho-C consists of two stars which are resolved 
only in HST images. It consists of a brighter bluer 
component {B = 21.28, V = 19.38, R = 18.10 
mag; RP04) and a dimmer redder component 
{B = 22.91, V = 20.53, R = 19.23 mag; RP04). 
In our analysis we find a consistent solution for the 
spectrum and infer that this is from the brighter 
bluer component. We find that Tycho-C is a 
metal-poor giant, probably located beyond the 
remnant. Tycho-C, similarly to Tycho-A, might 
be compatible with a giant-donor-star scenario. 
Its lack of rotation and kinematics, however, make 
it an uncompelling candidate. The only informa- 
tion we have about the dimmer component is the 
proper motion, which is insignificant with /x^ = 
0.58±1.73mas yr-\/X5 = -0.29± 1.21mas yr'^ 

Tycho-D is roughly a factor of 10 dimmer than 
the nearby star Tycho-C (separation « 0.6"). We 
could not measure reliable EWs for spectra with 
this S/N. Visual inspection of the star's spectral 
features shows it to be consistent with a cool star 
with low rotation. Its luminosity precludes it from 
being a relatively slowly rotating giant, and its 
slow rotation precludes it from being a siibgiant 
or main-sequence donor star. All of this suggests 
that Tycho-D is an uncompelling donor candidate. 

Tycho-E is the most distant star in this set 
(11.2 kpc), although large uncertainties in the dis- 
tance remain. It seems to be similar to Tycho-G in 
temperature, but appears to have a lower surface 
gravity. It is located 7" from the geometric center, 
but has no unusual stellar parameters or kinemat- 
ics. GII09 have suggested this to be a double- 
lined binary, but we are unable to confirm this 
using Fourier cross-correlation techniques. Ihara 
et al. (2007) have looked at iron absorption lines 
in stellar spectra made by the remnant and found 
Tycho-E to be unusual. They suggest that a star 
in the background would show blueshifted and red- 
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Fig. 6. — Observed spectra of Tycho-G centered around the Li A6707 line. Synthetic spectra with different Li 
abundances are overplotted. The thick red Une represents the Li abundance corresponding to the best-fitting 
value, and unsatisfactory fits (±0.15 dex) are plotted as thin black lines. 



shifted iron lines, whereas a star inside the rem- 
nant would only show blueshifted iron lines, and a 
foreground star would not show any iron features 
from the remnant. Ihara et al. (2007) claim that 
Tycho-E only shows blueshifted lines, and thus 
suggest that it is inside the remnant. We believe, 
however, that Tycho-E is located far behind the 
remnant and suggest that a low column density 
on the receding side of the remnant could cause 
a lack of redshifted iron features. In summary, 
a lack of rotation, kinematic signatures, and an 
inconsistent distance make Tycho-E a very weak 
candidate. 

Tycho-G is located 30" from the X-ray center, 
making it the most remote object from the center 
in this work (in the plane of the sky; for compari- 
son a distance of 32.6" corresponds to 1000 km s^^ 



over 433 yr at the distance of 2.8 kpc). This work 
confirms the radial velocity measured by GH09 
and WEK09. Figure 2 shows the expected dis- 
tribution of radial velocities from the Besangon 
model of Galactic dynamics. Tycho-G lies well 
within the expected range of radial velocity for 
stars with its stellar parameters and distance. 

In addition, this work has analyzed the proper 
motion of stars around the center of SN 1572. 
Figures 1 and 11 show Tycho-G not to be sig- 
nificantly deviant from the distribution of proper 
motions in the SNR 1572 neighborhood. Figure 1 
shows Tycho-G to be a 2a outlier, which implies 
that there should be about six stars in the HST 
sample sharing similar proper-motion features as 
Tycho-G; thus, its proper motion is by no means 
a unique trait. To further explore the proper- 
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Fig. 7. — EWs for the eight Ni lines in common between GH09 (open red squares) and this study (filled 
black circles) for Tycho-G. Lines (a-h) are 5082.35 A, 5088.54 A, 6086.28 A, 6175.37 A, 6176.82 A, 6643.64 A, 
7748.89 A, and 7797.59 A, respectively 



motion parameter space, we have selected candi- 
dates within a 1° radius around SN 1572 from the 
proper-motion catalogue PPMXL (Roeser et al. 
2010). To exclude the many foreground stars we 
introduced the additional selection criteria i? > 16 
mag and V — R < 1 mag (for comparison, the 
Sun has a color of V — R = 1.3 mag). These se- 
lection criteria are meant to exclude foreground 
stars. We tested this by applying the same se- 
lection criteria on the Besangon Model, resulting 
in 95% of stars more distant than 2kpc. This 
shows that a high proper motion at great distances 
is not a unique feature, as there are more stars 
that share this trait (see Figure 11). In particu- 
lar, stars in the thick disk have motions entirely 
consistent with Tycho-G (see contours in Figure 
11, and Figure 10 in GH09). Finally, the HST 
proper-motion measurements are challenging, and 



it is conceivable that there are systematic errors in 
our proper-motion measurements which are larger 
than our reported statistical errors. Such errors 
would tend to increase the chance of larger-than- 
actual proper-motion measurements. Taken in to- 
tal, while Tycho-G may have an unusual proper 
motion, the significance of this motion, even if 
current measurements are exactly correct, is not 
exceptional. 

As described, the kinematic features of a donor 
star might easily be lost in the kinematic noise 
of the Galaxy. WEK09 recommend using post- 
explosion stellar rotation as an additional possible 
feature for a donor star. This work suggests that 
Tycho-G has a rotation below the instrumental 
profile of 6 km s~^, much less than expected for a 
donor star (for an estimate, see Kerzendorf et al. 
2009). New results by Pan et al. (2012), how- 



26 



1.1 




6085 6086 6087 6088 6173 6174 6175 6176 

Wavelength(A) Wavelength(A) 

Fig. 8. — Observed spectra centered around five Ni lines in common with GH09 for Tycho-G. Synthetic 
spectra with different Ni abundances arc overplotted. The thick red line represents the Ni abundance 
corresponding to the value derived from EW analysis, and unsatisfactory fits (±0.3 dex) are plotted as thin 
black lines. 



ever, suggest that only taking tidal coupling into 
account could overestimate the rotation, and thus 
Tycho-G's low rotation might still be reconcilable 
with a donor-star model. 

We find Tycho-G to be a subgiant /main- 
sequence star with roughly solar temperature and 
metallicity. GH09 measure a nickel enhancement, 
which they believe to originate in the contam- 
ination from the ejecta. We have conducted a 
detailed comparison with GH09's measurement in 
§ 3.5 and do not find Tycho-G to be an outlier 
as suggested by GH09, but rather consistent with 
other stars of similar metallicity. In addition, our 



Li measurement is in agreement with that of GH09 
(see Tabic 6). In contrast to the GII09 interpre- 
tation, this Li abundance is consistent with that 
of stars of similar parameters (Baumann et al. 
2010). Finally, we have measured the distance 
to Tycho-G, showing it to be consistent with a 
background star. In addition, the radial-velocity 
signature matches that of background stars (see 
Figure 2). 

In summary, while Tycho-G may have unusual 
kinematics as indicated by its proper motion, the 
significance of this motion is not large when com- 
pared to a large sample of similar stars in the di- 
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Fig. 9. — Same as Figure 8 but for the remaining four Ni lines in common with GH09 (the upper-left line is 
also seen in the previous panel). 



rection of the Tycho remnant. Furthermore, such 
a kinematic signature, if it were related to the bi- 
nary orbital velocity, might predict rotation for 
Tycho-G which we do not observe (modulo the 
caveats from WEK09 & Pan et al. 2012). Also, 
we have not found a reasonable explanation for 
Tycho-G's large distance from the geometric cen- 
ter, and suggest that Tycho-G is unlikely to be 
related to the Tycho SNR. 

5. Conclusion 

This work did not detect an unambiguously 
identifiable donor-star candidate. Although Tycho-B 
shows some unusual features, there currently re- 
mains no convincing explanation for all of its pa- 



rameters which can be attributed to the donor- 
star scenario. We believe that our results provide 
evidence that the Tycho SNR does not have a 
main-sequence, subgiant, or red giant donor star. 
Some other possibilities remain. In the spin-down 
scenario, the companion star can become a helium 
white dwarf from a red giant donor, or a very low 
mass main-sequence star from a more massive 
main-sequence star. Such a compact compan- 
ion can escape detection (Di Stefano et al. 2011; 
Justham 2011; Hachisu et al. 2012a, b). Another 
scenario is a helium donor, such as the so-called 
sub-Chandrasekhar mass explosions discussed by 
Livne & Arnett (1995) and Sim et al. (2010). 
These progenitor systems might leave a very faint 
and fast-moving helium star, or no remnant at all 
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Fig. 10. — The figures show error contours for distance, extinction, and mass of the candidates. In the 
distance plots we indicate the distance range of SNR 1572 with a gray shade. The lower right shows the 
optimal isochrone (Pietrinferni et al. 2004) for the measured values of T^s and log g. 
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Table 10 

Distances. Ages, and Masses of Candidate Stars 



Tycho 


Mass 


CMass 


Age 


CTAge 


D 




(Name) 


(M/Mq) 


[M/Mq) 


(Gyr) 


(Gyr) 


(kpc) 


(kpc) 



Tycho-A 


2.4 


0.8 


0.7 


2.3 


1.4 


0.8 


Tycho-B 


1.8 


0.4 


0.8 


0.3 


1.8 


0.8 


Tycho-C 


0.9 


0.4 


10.0 


3.4 


5.5 


3.5 


Tycho-E 


1.7 


0.4 


1.4 


1.1 


11.2 


7.5 


Tycho-G 


1.1 


0.2 


5.7 


2.1 


3.7 


1.5 



(R. Pakmor 2012, privat communication). Such 
a progenitor would probably evade detection, and 
would likely not leave traces, such as circumstellar 
interaction with the remnant, or early hght-curve 
anomalies (Kasen 2010). Deep multi-epoch wide- 
field optical images should catch any such star 
speeding away from the remnant's center, but such 
observations have not yet been taken. Finally, a 
double-degenerate progenitor, in most cases, does 
not leave a compact remnant, and is consistent 
with our finding no donor star in SNR 1572. 

SN 1006 and SN 1604 (Kepler's SN) are two 
other SN la remnants in the Milky Way SN 1006 
is far from the Galactic plane and shows no signs 
of circumstellar interaction. Kerzendorf et al. 
(2012b) have studied this remnant and have not 
found any unusual star that can be explained with 
a donor-star scenario (consistent with this work). 
SNR 1604, while far from the Galactic plane, 
shows circumstellar interaction with its remnant, 
and has all the indications of what might be ex- 
pected from a single-degenerate scenario with an 
asymptotic giant branch donor (Chiotellis et al. 
2011). Observations of these remnants will better 
establish if there is a continued pattern to the un- 
usual stars in SN la remnant centers, or whether 
the lack of viable donor stars persists in multiple 
systems. 
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